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Rate  constants  and  product  branching  ratios  for  PO.*Clv~  ions  reacting  with  H  and  H?  were  measured  in  a 
selected  ion  flow  tube  (SIFT)  from  298  to  500  K.  P02C1",  P02C12",  POCl2“,  and  POCl3"  were  all  unreactive 
with  H2,  having  a  rate  constant  with  an  upper  limit  of  <5  x  10"12  cm3  s"1.  P02C12"  did  not  react  with  H 
atoms  either,  having  a  similar  rate  constant  limit  of  <5  x  10"12  cm3  s-1.  The  rate  constants  for  PO?Cl", 
POCl2",  and  POCl3"  reacting  with  H  showed  no  temperature  dependence  over  the  limited  range  of  298— 500 
K  and  were  approximately  10—20%  of  the  collision  rate  constant.  Cl  abstraction  by  H  to  form  HC1  was  the 
predominant  product  channel  for  P02Cr,  POCl2~,  and  POCl3",  with  a  small  amount  of  Cl  observed  from 
POCl2"  +  H.  Reactions  of  02  and  03  with  the  POC1"  products  ions  from  the  reaction  of  POC1?"  +  H  were 
observed  to  yield  predominantly  P03"  and  P02",  respectively.  POC1"  reacted  with  02  and  03  with  rate  constants 
of  8.9  ±  1.1  x  10"11  and  5.2  ±  3.3  x  10“10  cm3  s"1,  respectively.  No  associative  electron  detachment  in  the 
reactions  with  H  atoms  was  observed  with  any  of  the  reactant  ions;  however,  detachment  was  observed  with 
a  PO"  secondary  product  ion  at  high  H  atom  concentrations.  Results  of  new  G3  theoretical  calculations  of 
optimized  geometries  and  energies  for  the  products  observed  are  discussed. 


Introduction 

Recently,  reactions  of  various  positive  and  negative  ions  with 
H  and  H2  have  been  studied  to  ascertain  their  relevance  to  the 
space  environment.1"4  H2  and  H  are  important  species  in  com¬ 
bustion  processes  as  well.  Ion— molecule  chemistry  from  air 
plasmas  has  been  shown  to  be  a  viable  method  for  enhancing 
combustion.5’6  However,  the  role  of  negative  ion  chemistry  in 
combustion  is  not  as  well  understood.  In  addition,  oxyphos- 
phorus  compounds  have  been  shown  to  be  reasonable  alterna¬ 
tives  to  phosphine  additives  for  hydrogen-fueled  scramjet  com¬ 
bustors,  yielding  a  measurable  increase  of  the  flameholding  sta¬ 
bility  and  output  thrust  in  the  supersonic  expansion  because  of 
enhanced  H,  O,  and  OH  recombination  kinetics.7  Flames  doped 
with  oxyphosphorus  compounds  have  also  been  observed  to  pro¬ 
duce  anions  of  these  compounds  whose  profile  along  the  flame 
axis  is  strongly  influenced  by  reactions  with  H  as  well  as  H2.8 

The  chemistry  of  anions  generated  from  phosphorus  oxy¬ 
chloride,  POCl3,  has  recently  been  studied  in  our  laboratory 
using  a  selected  ion  flow  tube  (SIFT).9  POCl3  is  known  to 
generate  Cl",  POCl2“,  and  POCI3"  by  dissociative  electron 
attachment.10’11  Therefore,  the  rate  constants  and  product 
branching  ratios  of  PO*Clv"  ion  reactions  with  H2  and  H, 
including  P02C12"  and  P02C1"  ions  produced  by  source 
chemistry  from  P0C13,12  have  been  measured  in  the  SIFT  from 
298  to  500  K.  The  current  experiments  address  the  issues 
outlined  above  and  further  our  recent  studies  of  02  and  03 
reactions  with  these  ions.9 

Experimental  Section 

(a)  Kinetics  Measurements.  The  SIFT  instrument  at  the  Air 
Force  Research  Laboratory  has  been  described  in  detail  else¬ 
where13  and  will  only  briefly  be  described  with  regard  to  the 
current  experiments.  All  of  the  reagents  have  been  used  as 
obtained  from  the  manufacturer  except  for  POCl3,  which  was 
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subjected  to  several  freeze— pump— thaw  cycles  to  remove 
dissolved  gases.  Electron  impact  on  POCl3  (Aldrich,  99%)  in 
an  effusive  source  using  a  thoriated  iridium  filament  produced 
ions  which  were  mass  selected  by  a  quadrupole  mass  filter  and 
injected  into  a  fast  flow  of  helium  buffer  gas  (AGA,  99.997%). 
The  buffer  gas  was  introduced  into  the  flow  tube  through  a 
Venturi  inlet  and  entrained  the  injected  ions.  The  He  used  for 
the  buffer  flow  and  discharge  source  was  passed  through  a 
liquid-nitrogen-cooled  sieve  trap  to  remove  water  vapor  and 
other  trace  impurities  before  entering  both  the  flow  tube  and 
atom  source.  H  atoms  were  generated  by  a  microwave  discharge 
on  a  mixture  of  He  and  H2  (Matheson,  99.9999%)  in  a  Pyrex 
tube.  Typical  H2  dissociation  fractions  were  5—10%.  The 
reaction  with  H2  alone  was  studied  first  because  H  and  H2  were 
added  simultaneously  to  the  reaction  region.  The  remaining 
parent  ions  and  any  product  ions  were  sampled  through  a  blunt 
nose  cone  aperture  and  mass  analyzed  with  a  second  quadmpole 
mass  filter  and  detected. 

Product  branching  ratios  were  obtained  by  extrapolating  the 
measured  branching  ratios  to  zero  neutral  concentration  to 
minimize  the  effects  of  secondary  chemistry.  The  product 
branching  ratios  for  the  major  species  have  uncertainties  of 
±10%.14  The  H  atom  concentration  was  calibrated  by  the  known 
rate  constant  for  Cl"  reacting  with  H  of  9.6  x  10"10  cm3  s"1, 
which  does  not  vary  over  our  temperature  range.15’16  The  rate 
constants  for  the  H  atom  reactions  have  relative  uncertainties 
of  ±20%  and  absolute  uncertainties  of  ±30%,  slightly  higher 
than  for  stable  gases  such  as  H2,  which  have  relative  and 
absolute  uncertainties  of  ±15%  and  ±25%,  respectively.13  In 
addition,  the  ion  current  at  the  nose  cone  can  be  monitored  as 
a  function  of  both  the  H  atom  concentration  and  total  ion  counts 
to  determine  if  associative  detachment  occurred,  because  product 
electrons  diffuse  rapidly  to  the  walls  resulting  in  a  decrease  in 
nose  cone  current.17’18 

To  probe  for  reaction  pathways  to  P03"  and  P02"  ions 
previously  observed  in  POCl3  doped  flames,19  the  chemistry  of 


10.1021/jp045203e  CCC:  $30.25  ©  2005  American  Chemical  Society 


Published  on  Web  02/24/2005 


20060117  470 


V. 


2560  J.  Phys.  Chem.  A,  Vol.  109,  No.  11,  2005 


Midey  et  al. 


TABLE  1:  Rate  Constants  and  Branching  Ratios  for  the 
Reaction  of  PO^C^r  4-  H  and  Secondary  Product  Ion 
Chemistry11 


rate  constant, 

[A'c]  (xlO-9  cm3  s-1) 
branching  ratios 


reactants 

products 

^'298 K\ 

kJ  mol  1 

300  K 

500  K 

0.36  [1.9] 

0.37  [1.9] 

POjCl"  +  H  — 

P02"  +  HC1 

-222 

1.00 

1.00 

HOPOCi  +  e- 

7 

<0.005  [1.9] 

<0.005  [1.9] 

po2c.2-h— 

P02C1"  +  HC1 

-27 

No  Rxn. 

No  Rxn. 

0.16  [1.9] 

0.16  [1.9] 

POCT-  +  H  — 

pocr  +  HC1 

-45 

0.94 

0.96 

Cl"  +  HPOC1 

-6 

0.06 

0.04 

HOPCl2  +  e" 

-10 

0=PHC.l2  +  e- 

-3 

0.38  [1.9] 

0.37  [1.9] 

POCI3-  +  H  — 

POCl2“  +  HC1 

-275 

1.00 

1.00 

HOPCb  +  e- 

44 

secondary 

cheraistry 

0.089  [0.62] 

POCl-Oj  — 

PO3-  +  Cl 

-4.08 

>0.99 

P02"  +  CIO 

-116 

<0.01 

0.52  [0.92] 

pocr  +  o3  — 

P02"  +  C102 

-45 

-0.80 

POr  +  Cl  +  02 

-32 

POjCl-  +  02 

-240 

-0.20 

“The  measured  rate  constant  is  given  with  the  Su-Chesnavich 
collision  rate  constant  value  in  brackets,28,29  both  in  units  of  IQ-9  cm3 
s"1.  The  heat  of  reaction  at  298  K,  A H%9gK,  has  been  calculated  with 
G3  theory  using  the  values  of  Fernandez  et  al.9  and  those  given  in 
Table  2.  All  of  the  ions  below  are  unreactive  with  H2,  having  an  upper 
limit  to  the  rate  constant  of  <5  x  10-12  cm3  s-1.  See  text  for  details. 


POC1”  ions  that  were  observed  as  product  ions  in  the  current 
H  atom  experiments  has  been  studied  as  follows.  H  atoms  were 
added  at  an  upstream  inlet  to  produce  POC1'  from  POCl2-  + 
H.  First,  O2  separately  and  then  a  3%  O3  in  O2  mixture20  were 
introduced  at  a  downstream  inlet  to  react  with  the  POC1”  ions, 
and  the  rate  constants  and  relative  branching  ratios  were 
measured  for  O2  and  then  O3.  It  is  optimal  if  a  sufficient  amount 
of  H  atoms  can  be  added  to  drive  the  POCl2-  +  H  reaction  to 


completion  before  examining  the  secondaiy  chemistry  of  the 
products.  However,  the  reaction  conditions  could  not  be  adjusted 
to  achieve  this  condition  given  the  H  atom  concentrations 
achievable  with  the  current  source.  Consequently,  a  large 
POCl2-  signal  remained  along  with  a  large  Cl-  signal  arising 
from  the  POCl2-  +  O3  reaction,  where  Cl-  accounted  for  >75% 
of  the  reactivity.9  The  POC1-  signal  was  much  smaller  than 
both  the  remaining  POCl2~  and  the  Cl~  signal  just  described, 
meaning  that  any  increase  in  Cl“  would  have  been  a  small 
addition  to  a  large  peak  that  could  not  be  measured.  It  should 
also  be  noted  that  the  operating  conditions  of  the  mass 
spectrometer  needed  to  resolve  m/z  —  79  (PO3”)  vs  m/z  =  82 
(POC1-)  severely  discriminated  against  m/z  =  35.  Thus,  these 
conditions  precluded  determining  the  branching  ratio  of  Cl-  that 
might  be  generated  by  the  POC1-  reaction.  Nevertheless,  the 
mass  discrimination  against  m/z  =  63  (PO2  )  is  much  less 
severe,  and  a  correction  to  the  branching  ratios  for  the  02  and 
O3  reactions  was  subsequently  applied.  Thus,  the  branching 
ratios  for  this  secondary  chemistry  with  02  and  O3  neglect  any 
contribution  from  Cl-  products  from  these  reactions.  All  of  the 
experimental  results  are  given  in  Table  1. 

(b)  Structure  Calculations.  The  total  energy  at  0  K  and  the 
enthalpy  at  298  K  for  optimized  geometries  of  the  relevant 
molecular  species  have  been  calculated  with  G3  theory21,22  using 
Gaussian  03W23  with  the  results  shown  in  Table  2.  The  structural 
parameters  given  in  Tables  3  and  4  were  determined  from  the 
MP2(FuIl)/6-3  lG(c/)  optimization  step  ofG3.  Electron  affinities 
(EA)  were  also  determined  for  PO,  P02C1,  and  POC1  and  are 
compared  with  literature  values.9,10,24-26 

The  wave  function  stability  of  the  molecules  shown  in  Tables 
2-4  was  verified  for  the  Hartree-Fock  (HF)  level  optimization 
step  of  G3.  However,  two  exceptions  required  additional 
calculations  to  obtain  geometries  with  a  corresponding  stable 
wave  function.  A  stable  wave  function  at  the  HF  level  could 
not  easily  be  obtained  for  the  isomers  P— O— H  and  0=PH0C1, 
molecules  that  could  be  produced  via  associative  detachment 
of  H  to  PO-  and  P02C1“,  respectively.  Consequently,  G3(B3LYP) 
calculations27  were  performed  for  two  sets  of  related  species: 
the  POH  and  HPO  isomers  and  the  HP02C1,  0=P0HC1,  P02C1, 
P02C1-,  and  H  reactants  and  products.  Performing  the  additional 
calculations  at  the  same  level  for  all  of  the  related  molecules 
allowed  for  comparisons  of  the  relative  energies  and  enthalpies 


TABLE  2:  G3  Total  Energies,  Enthalpies,  and  Electron  Affinities  (EA)* 


total  energy  (0  K)  enthalpy  (298  K)  EA  (calc.)  EA  (lit.) 


po  (cs„,  2nr) 

-416.372258 

-416.368940 

1.103 

1.092,‘T.OO4 

PO-  (CB*  32-) 

-416.412774 

-416.409438 

0=POO  (C„  2A§ 

-566.653441 

-566.648080 

3.311 

0=POO“  (C„  >A§ 

-566.775100 

-566.769614 

POC1  (C„  ]A§ 

-876.499870 

-876.495344 

1.342 

1.06,“1.25l/ 

POCT  (C„  2A§ 

-876.549106 

-876.544227 

PO2CI  (C2„,  'A,) 

—951.716178“ 

-951.711027“ 

2.109“ 

2.1457 

P02cr  (Cs,  2A§ 

-951.793661“ 

-951.788053“ 

HP=0  (Cs,  'A§ 

-416.976138 

-416.972310 

HO— PCI  (Ci,  2A) 

-877.064143 

-877.059023 

0=PH— Cl  (Ci, 2 A) 

-877.050838 

-877.046188 

HO— PC12  (Cs, 

-1337.185569 

-1337.179430 

0=PH— CE  (Cs,  'A§ 

-1337.182583 

-1337.176990 

HOPClj  (C,, 2 A) 

-1797.214934 

-1797.206856 

HO-POC1  (C,, 2 A) 

-952.283854 

-952.278005 

-952.288108“ 

-952.282096“ 

0=PHOCl  (Cs,  2A§ 

-952.243667“ 

-952.238088“ 

H  (25,/2) 

-0.501003 

-0.498632 

-0.501087“ 

-0.498726“ 

HC1  (CB„,  >2+) 

-460.654664 

-460.651360 

“  Zittel  and  Linerberger.24  6  Wu  and  Tieman.25  “  Miller  et  al.,  MP4  theory.10  ‘'Miller  et  al.,  G2  theory.26  “  Calculated  using  G3(B3LYP)  method. 
See  text  for  details.  -''Fernandez  et  al.,  G3  theory.9  *  All  units  except  EA  are  in  hartrees.  EA  values  are  in  eV. 


Reactions  of  PO*Cl>’  Ions 


TABLE  3:  Bond  Lengths  (A)  and  Angles  (degrees)  for 
Geometries  from  G3  Theory 


parameter 

0=PC1 

0=PCF  HOPC1  0=PHC1  HP=0 

r  (P-O) 

1.493 

1.533  1.647  1.488 

1.517 

r  (P-Cl) 

2.078 

2.390  2.072  2.038 

r  (P— H) 

1.418 

1.453 

r  (O— H) 

0.976 

”  (Cl-P-O) 

110.5 

110.6  102.1  118.9 

”  (H— P— O) 

116.7 

105.6 

”  (H-P-Cl) 

99.4 

”  (P-O-H) 

114.7 

dihedral 

46.2 

”  (Cl-P-O-H) 

parameter 

HOPC1; 

,  HOPCf  0=PHC12  HOP=OCl 

r  (P—O1) 

1.618 

1.628  1.480 

1.489“ 

r  (P—O2) 

1.629 

r  (P— Cl1) 

2.194 

2.073  2.016 

2.041 

r  (P— Cl2) 

2.098 

2.073  2.016 

r  (P— Cl3) 

2.033 

r  (P— H) 

1.399 

r  (0— H) 

0.979 

0.979 

0.979 

”  (Cl'-P-O1) 

95.6 

101.4  115.7 

117.7 

”  (Cl'-P-O2) 

100.3 

”  (CF-P-O1) 

94.1 

101.4  115.7 

”  (CP-P-O1) 

105.3 

”  (Cl1— P— Cl2) 

147.6 

99.2  103.2 

”  (CP-P-CP) 

102.9 

’’  (CP-P-Cl3) 

104.2 

”  (H-P-O1) 

117.1 

’  (H-P-02> 

101.4 

”  (H-P-Cl2) 

101.4 

”  (P-O-H) 

111.5 

115.9 

109.9 

”  (O-P-O) 

116.0 

dihedral 

-19.3 

-51.0 

-143.8 

”  (Cl'-P—O— II) 

dihedral 

-168.3  51.0 

”  (CP-P-O-H) 

dihedral 

85.7 

”  (CP-P-O-H) 

dihedral 

-15.9 

”  (0'-P-02-H) 

“  P—O1  is  a  double  bond, 

P—O2  is  a  single  bond. 

TABLE  4:  Bond  Lengths  (A)  and  Angles  (degrees)  for 

Geometries  from  G3  Theory 

parameter 

P=0  P=0-  0=POO 

0=POO- 

r  (P—O1) 

1.538  1.595  1.487“ 

1.519“ 

r  (P—O2) 

1.704 

1.601 

r  (02-03) 

1.372 

1.514 

”  (O' -P-O2) 

105.2 

110.8 

"  (P-02-03) 

115.7 

98.7 

dihedral  ”  (O— P- 

O 

1 

O 

180.0 

180.0 

“  P—O1  is  the  terminal  P— O  double  bond. 


of  reaction  and  for  determination  of  stable  wave  functions,  where 
the  geometries  were  optimized  using  density  functional  theory 
(DFT)  at  the  B3LYP/6-3  lG(t/)  level  for  G3(B3LYP). 

The  wave  functions  of  these  oxyphosphorus  compounds  using 
G3(B3LYP)  were  all  stable  using  DFT,  with  the  only  exception 
still  being  the  high-lying  P— O— H  isomer.  However,  this  isomer 
lies  ca.  1.5  eV  higher  in  energy  than  H— P=0  according  to  both 
the  G3  and  DFT  calculations.  Although  the  energy  of  the  opti¬ 
mized  wave  function  is  thus  somewhat  uncertain,  the  P— O— H 
isomer  would  still  be  inaccessible  at  the  temperatures  studied 
here.  The  G3(B3LYP)  energies  are  also  given  in  Table  2. 

Results  and  Discussion 

(a)  H  and  H2  Reactions.  Rate  constants  and  branching  ratios 
for  the  reactions  of  P02Cr,  P02C12-,  POCl2”,  and  POCl3-  with 
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Figure  1.  Mass  spectra  for  the  reaction  of  POCF-  with  H  at  298  K  as 
measured  in  the  selected  ion  flow  tube  (SIFT).  The  dashed  line  is 
without  any  H  atoms  in  the  flow  tube  and  the  solid  line  is  with  H 
atoms  added  to  the  flow  tube. 

H  and  H2  measured  from  298  to  500  K  are  shown  in  Table  1. 
The  enthalpies  of  reaction  at  298  K,  A are  included  in 
Table  1  and  have  been  determined  from  both  the  previous9  and 
new  G3  calculations.  The  new  G3  calculations  of  the  total  ener¬ 
gy  at  0  K  and  the  enthalpy  at  298  K  are  presented  in  Table  2, 
representing  optimized  structures  given  in  Tables  3  and  4  as  men¬ 
tioned  above.  None  of  the  ions  currently  studied  react  with  H2 
to  a  measurable  extent;  therefore,  the  limitations  of  the  instrument 
permit  only  an  upper  bound  for  the  rate  constant  of  <5  x  10~12 
cm3  s-1  for  the  H2  reaction.  P02C12~  is  unreactive  with  H  atoms 
as  well.  A  comparable  limit  can  be  placed  on  the  rate  constant 
for  this  reaction  since  no  products  were  observed.  A  similar 
lack  of  reactivity  of  P02C12“  with  02,  O3,  and  even  POCI3  has 
been  previously  been  observed  in  the  SIFT.9'12  This  ion  is  very 
stable  as  evidenced  by  G3  calculations  (EA  P02C12  =  5.788 
eV)9  and  is  a  prominent  terminal  ionic  species  in  the  reaction 
of  POCl2~  and  POCF-  with  O3.9  The  rate  constants  for  the 
reactive  species  are  roughly  10—20%  of  the  Su-Chesnavichvich 
collision  rate  constant28’29  and  show  no  temperature  dependence 
over  the  temperature  range  of  298—500  K  presently  studied. 

By  far  the  most  prevalent  channel  for  all  three  reactive  ions, 
P02C1",  POCl3-,  and  POCl2~,  is  Cl  abstraction  by  H  to  form 
HC1.  Cl  abstraction  is  also  exothermic  for  the  P02C12-  reaction 
with  H  even  though  the  ion  proves  to  be  unreactive.  Figure  1 
shows  mass  spectra  for  the  POCl2-  reaction  taken  with  and 
without  H  atoms  added.  The  mass  spectrum  shows  almost 
exclusively  POCl2"  before  the  H  atoms  are  added.  At  high 
concentrations  of  H  atoms,  three  products  are  observed:  Cl-, 
POCI-,  and  PO-.  However,  PO-  is  exclusively  a  secondary 
product  formed  from  POCI-  reacting  with  H  as  shown  in  eq  1. 

POCF  +  H  —  PO-  +  HC1  4-  46  kJ  moF1 
—  CF  +  HPO  +  47  kJ  moF1 
—  HOPC1  +  e~  +  48  kJ  moF1  (1) 

It  is  not  possible  to  completely  deconvolute  the  secondary 
chemistry  because  the  reaction  of  POCF  with  H  appears  to  be 
near  collisional;  that  is,  the  secondary  chemistry  is  much  faster 
than  the  primary  chemistry.  It  is  clear,  however,  that  PO-  is 
formed  in  that  reaction,  whereas  CF  and  e-  may  also  be  formed 
as  well  as  shown  in  eq  1 .  The  reaction  of  POCl2-  with  H  also 
has  a  minor  primary  product  channel  giving  CF  and  HPOC1 
that  is  only  slightly  exothermic,  but  it  is  possible  that  the  CF 
product  signal  is  primarily  from  the  POCF  reaction  in  eq  1 .  In 
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addition,  PO-  has  been  observed  to  undergo  associative 
detachment  with  H,  making  it  hard  to  distinguish  if  POC1-  does 
as  well.  Nevertheless,  forming  POHC1  with  the  H  bound  to  the 
O  via  this  route  is  48  kJ  mol-1  (0.50  eV)  exothermic  based  on 
the  current  G3  theoretical  results  in  Table  2. 

An  interesting  observation  is  that  none  of  the  primary  PO.vCl;r 
ions  studied  undergo  associative  detachment  with  H  atoms.  The 
reaction  of  POCl2“  with  H  is  >3  kJ  mol-1  (0.03  eV)  exothermic 
for  either  HPOCh  isomer  as  shown  in  Table  1,  perhaps 
indicating  barriers  to  forming  the  neutral  product.  Previous 
studies  of  PO-  ion  reactions  have  found  that  PO-  was  very 
difficult  to  make  in  the  ion  source  of  the  SIFT  so  that  insufficient 
signal  was  available  for  measuring  the  rate  constant  for  the  H 
atom  reaction,  given  the  relatively  low  total  concentrations  of 
H  as  well.30  The  final  product  fractions  in  Table  1  have  been 
corrected  for  the  secondary  chemistry  of  the  product  ions. 

As  mentioned  above,  the  G3  results  in  Tables  2—4  show  that 
a  stable  minimum  energy  structure  for  POHC1  exists  having  an 
H  bound  to  the  O  atom  with  the  O— H  bond  slightly  displaced 
from  the  Cl— P—O  plane.  Another  stable  structure  of  HPOC1 
with  H  bound  to  the  P  atom  lies  0.36  eV  higher  in  energy.  All 
of  the  other  neutral  molecules  containing  H  shown  in  Tables  2 
and  3  also  have  the  lowest  energy  structure  with  H  bound  to  O 
(using  G3  theory).  The  only  exception  is  H— P=0,  which  is 
ca.  1.5  eV  more  stable  than  P— O— H  as  discussed  above.  For 
HPO2CI,  the  isomer  having  a  P— H  bond  is  1.21  eV  higher  in 
energy  using  G3(B3LYP).  Therefore,  this  isomer  will  not  play 
a  role  up  to  500  K.  Interestingly,  the  O— H  bond  that  forms  is 
~0.98  A  in  all  of  the  optimized  structures  shown. 

The  bond  lengths  and  angles  for  the  POC1  neutral  calculated 
here  agree  well  with  previous  experimental31  and  theoretical31-33 
values.  The  EA  for  POC1  of  1.34  eV  calculated  here  agrees 
better  with  the  G2  value  of  Miller  et  al.26  of  1.26  eV  than  with 
the  lower  level  MP2  value  of  1.06  eV.10  The  G3  EA  of  PO  of 
1.10  eV  is  in  good  agreement  with  the  experimental  value  of 
Zittel  and  Lineberger  of  1 .092  eV  from  photoelectron  spectro¬ 
scopy,24  values  which  are  both  higher  than  the  1.00  eV  value 
of  Wu  and  Tieman.25  G3(B3LYP)  calculations  for  PO2CI  and 
PO2CI-  give  an  EA  of  2.109  eV,  in  excellent  agreement  with 
the  G3  value  calculated  recently.9  The  present  bond  lengths  for 
PO  and  PO—  are  both  ca.  0.06  A  longer  than  the  literature 
values.34  The  molecular  parameters  for  HPO  are  in  good 
agreement  with  experimental35  and  higher  level  theoretical 
values,35’36  with  the  current  MP2  P-0  bond  length  being  ca. 
0.03  A  longer  than  these  values.35’36 

(b)  Pathways  to  P02-  and  PO3-.  The  flame  results  of 
Goodings  and  co-workers  show  that  PO3-  and  PO2-  are  the 
main  ions  observed  when  POCl3  is  added  to  an  atmospheric 
pressure  methane -oxygen  flame.19  One  potential  first  step  in 
the  process  is  electron  attachment  to  form  POCI3-  and  POCI2-. 
Since  the  flame  is  at  atmospheric  pressure,  the  former  may 
dominate.11  Our  combined  work  on  the  H,  02  and  O3  reactions 
shows  that  only  a  minor  route  exists  for  generating  PO3-  from 
the  reactions  of  POCl2-  and  POCI3-,  both  of  which  can  be 
produced  by  electron  attachment  to  POCI3  as  just  discussed.10’11 
For  example,  POCI3-  reacts  with  H  to  form  POCl2"  which  is 
also  formed  directly  in  electron  attachment.  However,  there  is 
only  a  minor  channel  in  the  reaction  of  POCl2-  with  ozone  that 
forms  PO3-  (ca.  1%). 

As  seen  in  Table  1,  POC1-  reacts  with  02  with  a  rate  constant 
of  8.9  ±  1.1  x  10-11  cm3  s-1  to  form  >99%  PO3-  at  298  K, 
with  a  trace  amount  of  P02-  also  observed.  Correcting  the  03 
branching  ratios  and  rate  constants  for  the  02  reaction  is  also 
required,  and  the  results  for  03  in  Table  1  reflect  these 


corrections.  After  making  this  correction,  as  well  as  correcting 
for  the  ~1%  PO3-  contribution  from  the  POCl2-  reaction  with 
O3,  the  branching  ratios  for  POC1-  reacting  with  O3  are  ca.  0.80 
P02-  and  0.20  P02Cr  at  298  K.  The  rate  constant  for  the 
POC1-  +  O3  reaction  is  5.2  ±  3.3  x  10-10  cm3  s-1.  This 
measurement  has  a  much  larger  error  limit  than  typical 
measurements  because  of  the  small  concentration  of  O3  and  the 
small  deviations  in  successive  rate  constant  measurements  that 
accumulate  when  correcting  for  the  02  reactivity.  Nevertheless, 
these  results  clearly  show  that  the  reaction  of  POCl2-  with  H 
is  the  first  step  in  a  pathway  to  both  P02-  and  PO3-. 

The  minimum  energy  structures  for  PO3  and  PO3-  have  been 
shown  to  be  planar  with  three  equal  P—O  double  bonds  in  the 
anion.9  New  G3  calculations  show  that  peroxide  forms  of  PO3 
and  PO3-  exist  having  stable  planar  z-shaped  structures  with 
the  parameters  given  in  Table  4.  This  neutral  PO3  isomer  is 
3.199  eV  higher  in  energy  than  the  ground-state  PO3,  whereas 
the  analogous  PO3-  isomer  is  5.044  eV  higher  than  the  planar 
ground-state  PO3-  anion.  Even  though  the  reaction  of  POC1- 
with  02  is  408  kJ  mol-1  (4.23  eV)  exothermic,  this  higher  energy 
form  of  PO3-  will  not  be  accessible  in  the  current  temperature 
range.  However,  this  species  may  become  important  at  higher 
temperatures. 

Conclusions 

The  rate  constants,  and  branching  ratios  for  the  reactions  of 
P02C1-,  P02C12-,  PC)C12-,  and  POCI3-  with  H  and  H2  have 
been  measured  in  the  SIFT  from  298  to  500  K.  None  of  the 
ions  react  with  H2/and  P02C12-  does  not  react  with  H  either, 
placing  a  limit  on  the  rate  constant  of  <5  x  10-12  cm3  s-1  for 
all  of  the  systems  where  no  reaction  was  observed.  The  rate 
constants  for  H  atom  reactions  are  <20%  of  the  collision  rate 
constant  and  do  not  vary  with  temperature  from  298  to  500  K. 
The  only  product  observed  is  Cl  abstraction  by  H,  except  for  a 
minor  reaction  pathway  with  POCl2-  that  produces  Cl-  and 
HOPC1.  Combining  the  present  results  with  the  O3  reactions 
studied  previously9  shows  that  a  pathway  from  electrons  to  P02~ 
and  PO3-  exists  for  methane— oxygen  flames  doped  with 
POCI3,19  where  the  H  atom  reaction  products  can  react  with  02 
and  O3.  In  the  future,  we  also  plan  to  examine  reactions  with  O 
atoms  and  02(a  1  Ag)  for  more  robust  pathways  to  PO3-  and 
P02-. 
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